Mon. Not. R. Astron. Soc. 000, [[HH] (2011) Printed 9 January 2012 (MN lAT^X style file v2.2) 

SZ power spectrum and cluster numbers from an extended 
merger-tree model 



o 

(N 

in 



O 

U 

6 



> 
o 



o 



X 



Irina Dvorkin 1 *, Yoel Rephaeli 1 ' 2 , Meir Shimon 1 

1 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 

2 Center for Astrophysics and Space Sciences, University of California, San Diego, La Jolla, CA 92093-0424 



9 January 2012 



ABSTRACT 

Wc have recently developed an extended merger-tree model that efficiently follows hi- 
erarchical evolution of galaxy clusters and provides a quantitative description of both 
their dark matter and gas properties. We employed this diagnostic tool to calculate 
the thermal SZ power spectrum and cluster number counts, accounting explicitly for 
uncertainties in the relevant statistical and intrinsic cluster properties, such as the halo 
mass function and the gas equation of state. Results of these calculations are com- 
pared with those obtained from a direct analytic treatment and from hydrodynamical 
simulations. We show that under certain assumptions on the gas mass fraction our 
results are consistent with the latest SPT measurement. Our approach can be partic- 
ularly useful in predicting cluster number counts and their dependence on cluster and 
cosmological parameters. 

Key words: galaxies: clusters: general - large-scale structure of Universe - cosmic 
background radiation 



1 INTRODUCTION 

The importance of the Sunyaev-Zel'dovich (SZ) effect 
as a valuable cosmological probe and a usefull tool for 
discovery of high-redshift clusters is now firmly estab- 
lished. In view of the sensitivity of the effect to the 
, halo mass function (MF), in particular to the normal- 
ization of the matter power spectrum, the SZ effect can 
be used to constrain various cos mological parameters of 
the ( s tandard) ACDM model (e.g., Holder. Haiman fc Mohri 

~~ir 



l200ll : lAllen. Evrard fc MantdboTH . and references therein), 



including the eq uation of state of dark energy (e.g.. 



Wan g et al. 12004) and the to tal neutrino mass (e.g.. 



Shimon. Sadeh fc R ephaeli 2011), as well as alternative cos- 



mological models (e.g., Sadeh. Rephaeli fc S ilk 2007). 

However, the SZ effect is also sensitive to the intraclus- 
ter (IC) gas physics, which is poorly known at high redshifts. 
The latest observatio nal results from th e Atacama Cosmol- 
ogy Telescope (AC T: iFowler et al.ll20ld) and the South Pole 
Telescope (SPT: iLueker et all l20ld ; IShirokoff et~afl l201ll : 



iReichardt et al.ll2011 ) underline the need for realistic model 
ing of IC gas properties as a prerequisite for use of clusters as 
probes to determin e cosmological pa r amete rs from SZ mea- 
surements. Indeed, IReichardt et al.l |201 ll ) concluded that 
the theoretical uncertainty in predicting the SZ power spec- 



trum is significantl y larger than the statistical errors. Several 
theoretical studies JShaw et al.ll20ld : lTrac. Bode fc Ostrikerl 
l201ll : lEfstathiou fc Migliaccidl2oill ) have also demonstrated 
the sensitivity of the SZ effect to simplified modeling of IC 
gas, particularly its equation of state. 

There are currently two approaches to modeling the SZ 
powe r spectrum: An analytical approach (|Komatsu fc Seliakl 
2002) assumes universal profiles of the dark matter (DM) 
density and the IC gas density (and pressure), whose pa- 
rameters follow simple scaling laws with mass and the clus- 
ter redshift (of observation). The relative contribution of 
clusters at different redshifts is weighted according to an ap- 
propriate MF. The computational efficiency of this method 
makes it feasible to explore a large portion of the parame- 
ter space, but it relies heavily on average scaling relations 
between the cluster observables, such as the concentration 
parameter, gas temperature, and the cluster mass. In addi- 
tion, this description does not include the intrinsic scatter 
in these scaling relations which may depend on the cluster 
mass and redshift. 

The second approach i s based on numerical simula- 
tions , either dynamical (e.g.. IShaw et _a l. 2009; Schg al et al 
2010f) or hydrodynamical (e.g. , ISeliak. Burwell fc Pen 



200ll: ISpringel White fc Hernquistll200ll: iBond et all 12005 



E-mail: irina@wisc.tau.ac.il 



Schafer et al.l l200d : iRoncarelli et all 120071 : iBattaglia et all 
20ld ). from which large catalogues of clusters are created. 
This method reproduces the great variety in the observable 
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properties of clusters, but the large computational costs lead 
to two shortcomings: the volume of the simulation is lim- 
ited, which leads to insufficient statistical sampling of high- 
mass (~ 10 15 Mq) clusters, and it is highly inefficient in 
testing various cosmological and IC gas models. A varia- 
tion on the numerical approach uses the Lagrangian per- 
turbation theory (of which the Zel'dovich approximation is 
the first term) to generat e DM halo catalogues with much 
less computational effort jMonaco. Theuns fc Taffonill2002l ; 



iHolder. McCarthy fc Babulll2007l ). However, this technique 
too is limited by the size of t he simulation box. 

In this work we use our (|Dvorkin fc Rep haeli 20 ll]) ex- 
tended merger-tree model of cluster evolution to calculate 
the thermal SZ power spectrum and cluster number counts. 
This method enables quantitative predictions of the outcome 
of numerical simulations by taking into account the intrin- 
sic scatter in cluster properties that result from their differ- 
ent formation histories. Our model is an improvement over 
standard analytical methods in that it accounts for the full 
merger history and formation redshift of each cluster, and 
does not depend on pre-calibrated scaling relations. On the 
other hand, our approach can be used in studies of cosmolog- 
ical parameter estimation in a much more computationally 
efficient way than numerical simulations. 

This paper is organized as follows. In Section [5] we re- 
view our merger-tree model of cluster evolution and describe 
how it is integrated into the calculations of the SZ power 
spectrum and number counts. In Section[3]we present our re- 
sults and discuss their implications in Section [4] We use the 
following cosmological parameters: fi m — 0.25, Qa = 0.75, 
tt b = 0.045, H = 73 km/s/Mpc, a 8 = 0.8. 



2 EXTENDED MERGER-TREE MODEL 

We build merger tr ees of DM halos using the m odified 
GALFORM algorithm (|Parkinson. Cole fc HellyJ 200Sh which 
is ba sed on the excursion set formalism ( Lacev fc Cold 
1 19931 ). The conditional MF used in this algorithm is cal- 
ibrated to the outc ome of the Millennium Simulation 
jSpringel et al.l 120051) and is consistent with the Sheth- 
Tormen MF (|Sheth fc Tormenlll999r i. Consequently, we use 
this MF i n all of our ca l culati ons. The more updated MF pre- 
sented in iTinker et all (|2008T ) differs from the Sheth-Tormen 
MF by up to ~ 20% at z = 0, and by a higher fraction at 
larger redshifts. The adoption of a specific (especially an 
analytic, theoretically-based) MF, introduces modeling un- 
certainties; these are discussed in Section [3] 

For a cluster with a given mass and at a given observa- 
tion redshift we build a merger tree which represents its pos- 
sible evolutionary track. We define the major merger events 
as those with a mass ratio M> /M< < q for some q whose 
value is to be determined, and treat all other processes of 
mass growth as continuous (relatively) slow accretion. We 
assume that during the violent major merger events DM is 
redistributed in the cluster potential well, while the minor 
mergers and accretion processes affect mainly the outskirts 
of the cluster but not its central spherical region. We assume 
that the scale radius r s of the halo remains constant during 
slow accretion of matter onto the cluster, so that the radius 
of the halo R grows during accretion but the interior region 
remains essentially unchanged. Note that we deviate here 



some what from our original treatment (I Dvorkin fc Rephaelil 
1201 ll ) in which we assumed that the concentration parame- 
ter, c, is constant during slow accretion. 

Having built the merger tree, we start at the highest 
redshift with the smallest masses and calculate the density 
profile of each halo in the tree. To conform with common 
practice, and for direct compariso n with previous work, we 
assum e a Navarro-Frenk- White (jNavarro. Frenk fc White] 
1995) profile for all halos at all times: 



Pd{x) 



4p s 



x(l +x) 



(1) 



where x = r/r s , with r s the characteristic scale radius of the 
profile, and p s ~ p(l). Starting with an initial distribution 
of concentration parameters c(M, z) for the earliest halos, 
we calculate the concentration parameter of each successive 
halo in the tree from considerations of energy conservation. 
The outcome of this calculation is the concentration param- 
eter of the halo we started with at the specified redshift of 
observation. Generating a large number of trees gives an es- 
timate of the probability density function (PDF) for a given 
mass and redshi ft. Further details on the m erger-tree model 
can be found in iDvorkin fc Rephaelil l|201ll) . 

The next step is to model IC gas, which is assumed 
to constitute a small fraction of the cluster mass, and to 
not significantly affect the evolution of the cluster. The gas 
equation of state is assumed to be polytropic, namely that 
the pressure and density are related by 



Po(p/po) 



(2) 



with F = 1.2. 

The solution of the equation of hydrostatic equilibrium 
for a polytropic gas inside a p otential well of a DM halo is 
l|Ostriker. Bode fc Babu]|2005l ): 



p(x) 



po 



1 



B 



1 + n 



1 - 



ln(l + a;) 



where n — (T — 1) 1 , B is given by: 



B 



A-KGp s r s pm p 
k^To ' 



(3) 



(4) 



and fiirip is the mean molecular weight. The temperature 
profile is given by: 



T(x) = To 



1 



B 



1 - 



ln(l + x) 



(5) 



As a boundary condition we assume that the gas pres- 
sure at the virial radius obeys P g — f g Pd, where f g is the gas 
mass fraction and Pa = pd<T 2 . We obtain a 2 , the DM (3D) 
velocity dispersion, by solving the Jeans equation (in the 
field of NFW-distributed DM) . For the gas mass fraction we 
adopt a (n ominal) mean value of f g =0.1, estimated from 
the data in lBonamente et al.l (|2008T l. although we note that 
these measurements refer to the mass inside -R2500, within 
which the mean density is 2500 the background density (at 
the cluster redshift). While this may be considered a rela- 
tively small inner (and perhaps unrepresentative) region of 
the cluster, a similar mean va lue can be deduced from X-ray 
measurements within -R500 bv lVikhlinin et al.l (|2009T ). Below 
we explore other assumptions about f g . 

We also consider an alternative model for the IC gas, 
the fi— model with j3 = 2/3, in which case the density is 
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p{x) = 



Po 
1 + x 2 



(6) 



and the temperature is given by the solution of the equation 
of hydrostatic equilibrium with the (approximate) bound- 
ary condition T — > at large radii. Note that the non- 
dimensional radial coordinate is generally different from that 
in eq. flTJl. 

The SZ power spectrum is compu ted using the halo ap- 
proximation ()Komatsu fc Seliakll2002l ): 



c e = s( x y 



d ^ldz " 

dz 



dM^\y e (M,z)\ 2 



(7) 

where s(x) is the spectral dependence of the SZ signal 
given by: 



s(x) = X 



e x + 1 
e x - 1 



(8) 



where \ = hu/ksTo, V(z) is the comoving volume per stera- 
dian, dn/dM is the MF, and ye is the 2D Fourier transform 
of the projected Comptonization parameter, 



dXX tx/l. C(Z) 



(9) 



where c — R v /r 3 is the concentration parameter (R v is the 
virial radius), l s = dA(z)/r s , dji^z) is the angular diameter 
distance to the cluster, and ((x) is the gas (normalized) 
pressure 



COO 



^n e {x)T e {x). 



(10) 



Typical parameters are z max — 2, M min = 10 13 h 1 Mq and 
M max = 10 16 h^Me. 

We build one tree for each mass and redshift in equation 
(Q), and sum over all the mass and redshift range to obtain 
the power spectrum. When we repeat this calculation, the 
concentration parameter of each halo is slightly different, 
since it is effectively drawn from a probability distribution, 
resulting in a different power spectrum. However, for a suf- 
ficiently large number of mass and redshift bins that are 
used to calculate the integral in equation (Q), the calculation 
converges. Below we present the results of N tree s = 5 subse- 
quent calculations for each I. It is remarkable that the calcu- 
lation converges for such a small Agrees in comparison with 
the number of trees that were needed to estimate the con - 
centration parameter PDF in iDvorkin Sz Rephaelil |201l| y 
This is due to the fact that the concentration parameter 
PDF changes slowly as a function of mass and redshift, so 
for a dense grid of M and z in the evaluation of equation 
we adequately sample each PDF. 

We have also calculated th e expected cluster numbe r 
counts, following the method of ISadeh fc Rephaelil l|2004h . 
The SZ flux of a cluster is given by: 



af v = / R s (n,a B )Ai v (n)dn, 



(ii) 



where Cl is the direction on the sky, R s describes the detector 
beam with beam size ob and AI V is the (spectral) intensity 
change, 



A/„ = 



2(fc B T ) 3 
(he) 2 



■ y ■ g(x)- 



(12) 



The spectral dependence is g(x) = X e x (e x — 1)~ s(\), 
and To is the CMB temperature. Then the expected number 
of clusters with flux greater than some threshold AF„ is: 



N(AF V ) 



f ^Ldz f B(M,z) 
J dz J A p 



dN 
~dM^ 



dM. 



(13) 



The selection function is B(M,z) = 1 if the flux of a 
given cluster is above the threshold; otherwise, B(M, z) — 0. 
The calculations were carried out with Planck HFI 143 GHz 
channel, with a beam size of ob = 7.1', Av/v = 0.33, and 
flux sensitivity threshold of AF U = 12.6 mjy for 14 month 
observation. 

The computation time obviously depends on the num- 
ber of trees we grow for each mass and redshift, and the 
number of mass and redshift bins used to approximate the 
integral in equation (0. Calculating the power spectrum for 
different numbers of trees, and with various mass and red- 
shift resolutions, we find that the calculation converges for 
Ntmes = 5, A^AAf = 100 and Naz = 100, in which case it 
takes ~ 3 minutes to calculate the power spectrum Ce for a 
given I. Below we present results for 18 values of t between 
I = 100 and i = 9000. 



3 RESULTS 

The SZ power spectrum at v = 153 GHz, computed using 
our merger tree model, is shown in Figure Q] For comparison 
we show results of a standard calculation where we use the 
scaring relation fo r the concentration parameter deduced by 
I Duffy et al.l ()2008h from a sample of relaxed halos. Also plot- 
ted is the result of the standard calculation where we used 
the average c(M, z) relation from the merger-tree model: 

M 



log c(M,z) = 2.417 - 0.1081 • log ( — 



M 



- 1.57 ■ log(l + z) + 0.1039 ■ log I— J log(l + z) 

- 0.2486 -log 2 (l + z). 

(14) 

The normalization of the power spectrum is known to 
be very sensitive to as; we find the scaling Asz oc cr|' 4 , in 
accord with previous studies. The grey bands in Figure [T] 
correspond to erg = 0.8 ± 0.03, a level of uncertainty deter- 
mined from a joint an alysis of WMAP7, BAO and SN data 
l|Komatsu et al.ll201lh . 

The fiducial model presented in Figure [1] results in an 
overall normalization well above the level obtaine d from the 
latest SPT measurement dReichardt et al.ll201ll ): A S z(t = 
3000) = 3.65 ± 0.69 /iK 2 . Below we will show that within 
the uncertainties in the gas model, the gas mass fraction 
and the value of as, our calculation is consistent with these 
measurements. 

To compare the above results with those from simu- 
lations, we show in Figure [1] the power spectrum deter- 
mined from an adiabatic hydrodynamical simulation by 
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Figure 1. SZ power spectrum at v = 153 GHz. Results shown are 
from the standard calculation (thin blue line), merger-tree calcu- 
lation (thick pink line), a standard calculation with c(M,z) fit 
to the merger-tree model ( dot-dashed red line), ad iabatic hydro- 
dynamical simulations from Battaglia et al. ( 2010:) (dashed black 
line) . The grey area shows the uncertainty of the merger-tree cal- 
culation due to variations of ag. 



iBattaglia et al.l |2010T ) Q rescaled to v = 153 GHz. It is 
very interesting that both the normalization and the shape 
of the power spectrum obtained in our treatment agree well 
with the full numerical calculation, in spite of the very differ- 
ent approach used in the simulations. IBattaglia et al.l (|2010l ) 
also calculated the power spectrum in the case when ra- 
diative cooling, star formation, supernova and AGN feed- 
back were included, and found a reduction of power at 
small scales due to expansion of the gas to l arger radii. 
Similar resu l ts wer e obtained by IShaw et al.1 |2009l ) and 
ISehgal et al] l|20lut ). with different normalizations depend- 
ing on the amount of energy feedback. While we have not 
attempted to model feedback mechanisms in this work, the 
agreement of our model with numerical calculations in the 
simple adiabatic case shows that our description of the hier- 
archical growth of clusters and the intrinsic scatter in their 
properties is reasonably realistic. We stress that even though 
we make several simplifying assumptions, including that of 
hydrostatic equilibrium, we still get a better agreement with 
simulations than in the standard approach. The modeling of 
IC gas can be further refined, as we discuss below. 

There are two major differences between the merger- 
tree calculation a nd the standard tr eatment with the scal- 
ing relation from IDuffv et all l|2008h : the normalization in 
the former case is higher, and the peak shifts to higher mul- 
tipoles (£ ~ 6000 vs. i ~ 4000 in the latter case). These 
features can be understood in terms of the modified c(M, z) 
relation that is deduced from the merger-tree algorithm. 
When we fit the average c(M, z) relation from the merger- 
tree model to a general form c = A(M/M P i VO t) B (l+z) c that 
is used bv lDuffv etahl i|2008r i, we obtain A = 11.99tg:«, 
B = -0.079 ± 0.005 and C = -0.34 ± 0.03 for M pivot = 
2 x IO^/i^Mq, with 2a confidence level (CL), but this 
function provi des a poorer fit to our results than equation 
(|14[) . whereas IDuffv et al.1 (|200St ) obtain for their sample 



1 http:/ /www. astro. utoronto.ca/~battaglia/ 



B = -0.09 ± 0.009 and 
C = -0.69 ± 0.05 (la CL). Thus, the merger-tree model 
predicts a slightly higher normalization of the c(M, z) re- 
lation, which is partly responsible for the higher normal- 
ization of the SZ power spectrum, and a significantly lower 
redshift dependence, which also contributes to the increase 
in the power spectrum normalization and shifts the peak 
to higher multipoles (since clusters at higher redshifts are 
more concentrated and have a larger contribution to the 
SZ power). Indeed, when we calcul ate the SZ power spec- 
trum with the scaling relation from lDuffv et all (|2008h . but 
with the redshift dependence parameter C = —0.34 taken 
from our merger-tree model, the peak shifts to I ~ 5000 and 
the normalization increases by ~ 14%. An additional dif- 
ference between the calculation performed with the c(M, z) 
relation from equation (|14p and the full merger-tree calcu- 
lation, appears because the concentration parameter has an 
approximately log-normal distribution function: there are 
more halos with higher-than-average c than with lower-than- 
average. 

To what extent are these predictions for the concen- 
tration parameter reliable? There is some observational ev- 
idence that the concentration parameter is in fact higher 
than predicted by N-body simulatio ns. For example, the con - 
centration parameters measured bv lSchmidt fc Alien! (|2007T ) 
from X-ray observations of 34 dynamically relaxed clus- 
ters in the redshift range z = 0.06 — 0.7 are sign ificantly 
greater than those predicted by IDuffv et all i|2008l ) for the 
same masses and redshi ft s. In a large dataset compiled by 
IComerford fc Natara ian (2007), which consists of 62 clus- 
ters at redshifts up to z = 0.89 and masses up to at least 
~ 2 x 1O 15 M0, the normalization of the concentration pa- 
rameter is found to be higher by at least 20% than the results 
of numerical simulation s. A similar conclusion is reached by 
IWoitak fc Lokasl (|201Ch . who analyzed kinematic data of 41 
clusters at redshifts z < 0.1. A recent X-ray analysi s of 44 
clusters in the redshift range 0.1-0.3 bv lEtt'ori et aLI (|2010h 
found a normalization of c consistent with numerical simu- 
lations, although note that their constraint on the cosmo- 
logical parameters produces a rather high erg = 1.0 ± 0.2. 

Interestin gly, several strong lensing studies of massive 
clusters (e.g., IZitrin et al.l 120111 ) found anomalously large 
Einstein radii and concentration parameters when compared 
with predictions from numerical simulations. The results of 
the merger-tree model suggest that perhaps these concentra- 
tion parameters, although higher than the average values, 
are not incompatible with the ACDM model. 

It is important to note that the simulations of 

li 



im port a i 

IDuffv et al.l ([2008) span the mass range of 10 
10 15 h~ x Mq, similar to other numerical studies, while the 
masses that are important in the context of the SZ effect 
arc 10 13 -10 16 ft- 1 M Q . Low -mass and low-redshift halos out- 
number high-mass halos by several orders of magnitude; the 
c(M, z) relations are thus heavily weighted by the statistics 

of the low mass halos. 

For example, the sample of lDuffv et al.l |2008h contains 
1269 sufficiently well-resolved halos at z — 0, of which typ- 
ically less than ~ 0.1% are above 1O 14 /i _1 M0 if their rel- 
ative abundance approximately follows the Sheth-Tormen 
MF. At z = 0.5 only ~ 0.03% of all the halos are expected 
to have masses above lO 14 fe _1 M0. On the other hand, 
the parameters B and C are expected to vary with both 
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Figure 2. SZ power spectrum for different density profiles of DM 
halos. Shown are results for an NFW profile (solid line) and the 
profile in equation 1151 with a = 3 (dashed line). 
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Figure 4. SZ power spectrum with merger-tree model parameters 
in the range q = 7 — f3, k = 4 — 6 (grey area) and the fiducial 
model q = 10, n = 5 (thick line). 
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Figure 3. SZ power spectrum for different gas models: polytropic 
with different values of the adiabatic index, non-thermal pressure 
component, and /3-model with /3 = 2/3. 

mass and redshift llGao et al 1 120081 : |m 

unoz -Cuartas et all 

l201ll ; IPrada et al.ll201lh . Indeed, when we fit the average 
merger-tree results only in the range z = — 1, M = 
10 13 - lO 14 ^" 1 ^ we obtain B = - 0.09, C = -0.46, closer 
to the values in iDuffv et al.l |2008t ). This example demon- 
strates the inherent difficulty of numerical simulations to 
produce large samples of halos due to computational lim- 
itations. It is possible that because of the relatively small 
sample size numerical simulations overestimate the redshift 
evolution of the concentration parameter which biases the 
predicted SZ power spectrum. 

Our merger-tree approach enables us to test alternative 
DM halo density profiles, which is not possible in analytical 
calculations that are based on pre-calibrated scaling rela- 
tions. As an example, we explore the following modification 
of the NFW profile: 

Pdir) = r/r.(lVr/r.) a (15) 
for a > 2 (equations (|3][5]) are replaced by the appropriate 



solution for the new DM potential). Figure [2] shows that 
this modification with a = 3 has a minor effect on the am- 
plitude of the power spectrum, as it affects only weakly the 
integrated gas pressure in the outer regions of clusters. 

All the calculations in this work were performed with 
the simple IC gas model described i n Section [2] As 
was recently shown in several studies JShaw et al.(l201ol : 
iTrac, Bode fc Ostrikerll201ll : lEfstathiou fc MigliaccidboilT ). 
the strength of the SZ effect is very sensitive to the details of 
IC gas physics, for example a modest amount of non-thermal 
pressure can significantly lower the power spectrum. Here we 
explore the possibility of non-thermal pressure support, as 
well as slight variations in the parameters of the polytropic 
model. Figure [3] compares the results for different values of 
the adiabatic index T, a /3— model and a polytropic model 
with a non-thermal pressure component. For simplicity we 
assume that the non-thermal component is a constant frac- 
tion of the total pressure: 

Ptot — Pthermal + Prel — (1 + Srel) Pthermal , (1°) 

with S re i —0.1 as a representative value. 

It can be seen that the SZ power spectrum is very sen- 
sitive to slight changes in the value of the adiabatic index 
away from V — 1.2, an average value deduced from numeri- 
cal simulations. In accord with previous studies, we find that 
a small constant fraction of non-thermal pressure obviously 
reduces the normalization of the power spectrum but does 
not affect its shape. 

Figure U shows the uncertainty in our calculation due 
to the merger-tree model parameters q, which is the max- 
imal allowed ratio for major merger events, and k, which 
describes the initial distance between two clusters that are 
about to merge. As expected, this uncertainty is smaller 
than ~ 7%. 

Our fiducial model assumes f g = 0.1, a constant gas 
mass fraction across the whole mass and redshift range. 
However, the fraction of hot gas in clusters varies with mass 
and redshift of observation. The mass dependence arises 
from various galactic processes including star formation, 
ram pressure stripping of gas from galaxies, and galactic 
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Figure 5. SZ power spectrum calculated with the gas mass frac- 
tion taken from equation l|17|l and 5 re i = 0,1. Blue dot with er ror 
bars represents the SPT measurement iReichardt et al1i201ll) at 
I = 3000. The grey area shows the uncertainty due to variation 
in cr 8 . 



Figure 6. SZ cluster number counts above a specified redshift. 
Shown arc results of a s tandard calculation with a c[M, z) relation 
from lDuffv et al.l ll2008f ) (red dot-dashed line), with c(M, z) from a 
fit to the merger-tree model (blue solid line), and the full merger- 
tree calculation (black lines). 



winds. While the full redshift evolution of the gas mass 
fraction is not known, it clearly reflects aspects of cluster 
formation and evolution, which include also the effects of 
internal processes that re-distribute cluster baryons, those 
same processes that also imprint the (related) mass depen- 
dence. The latter was deduced in several X-ray studies of low 
redshift clusters JVikhlinin et alJIiobH : iGiodini et ai1l2009[ ) . 
iGiodini et al. (2009) derived the following dependence based 
on a combined sample o f 41 clusters at z < 0.2 observed 
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(17) 

We note that this fit is not valid beyond M ~ 1O 15 M . 

Not knowing the explicit scaling of the gas mass frac- 
tion with redshift, we can only account for its mass depen- 
dence by adopting the relation f g = M gaa , 50 o/Mtot,500 = 
M gaStViria i/Mtot,virial- To avoid having a baryon fraction 
greater than the cosmic mean, which arises in high-mass 
clusters if the scaling relations above are strictly followed, 
we set an upper limit of f cosmic = flb/Qm- 

The power spectrum, calculated with the gas mass frac- 
tion from equation (|17[) and assuming 8 re i = 0.1 is shown 
in Figure [5] It can be seen that our result is consistent with 
the SPT measuremen t if we allow for an uncertainty in <rg 
l|Komatsu et al.ll201ll ): as = 0.8 ± 0.03. 

While we did not account for the possible redshift evo- 
lution of the gas mass fraction, we note that precise mea- 
surements of the SZ power spectrum may actually be used 
to deduce, or significantly constrain the redshift evolution 
of f g , especially so when as is more precisely determined. 

The expected cluster number counts are presented in 
Figure [6] We show t he results of a stand ard calculation with 
c(M, z) taken from iDuffv et al.l {2008) and from equation 
(|14p . as well as the merger-tree calculation repeated 10 times 
so as to probe the distribution functions of the cluster prop- 
erties. There is ~ 15% difference between the results of the 



merger-tree model and the standard calculation, a level of 
variance that could have appreciable ramifications on precise 
cosmologi cal parameter extraction based on cluster numb er 
counts (|Holder. Haiman fc Mohrll200ll ; IWang et al.ll2004l ). 

A major source of uncertainty in calculations of the SZ 
effect is the halo mass function. It is evident that the ac- 
curacy with which the MF can be determined in a given 
numerical simulation is reduced at higher redshifts, espe- 
cially at the high mass end, due to the limited volume of 
the simulation. To illustrate this point, we can estimate the 
expected number of halos at different redshifts according to 
the Tinker MF in the simulations that were used to calibrate 
it. For example, for z = 1.25 the number of halos per unit 
In M is expected to be less than ~ 100 for M > 3 ■ 10 14 M Q , 
while at z = 2.5 there is less than one halo with mass around 
2-lO 14 M0 per unit In M in the combined simulation volume. 

It can be argued that large masses are rare at high red- 
shifts and their contribution to the SZ power spectrum, al- 
though uncertain, is not important. Figure [7] shows the rel- 
ative contribution of different masses and redshifts to the 
integral in equation Q in the following way: we calculate 
this expression for £ = 4000 with different z max (x-axis) and 
Mmax (y-axis) and plot the result as a fraction of the fiducial 
value for which we took z max = 2 and M max — 10 16 /i _1 Mq. 
The contours are for constant fraction of the fiducial valuqfl- 
Two important features can be seen in Figure [7] (a) we have 
to integrate at least up to z ~ 1.4 to have a reliable estimate 
of the effect, and (b) masses above 2 • 10 14 Mq contribute 
roughly ~ 50% of the effect (over the whole redshift range) . 
It is clear that a large error in the MF even at z ~ 1.25 
(which only worsens at higher redshifts) is very detrimental 
to the reliability of the calculation. 

Finally, it should be emphasized that there is an inher- 
ent limitation in the accuracy of cluster MF deduced from 
dynamical cosmological simulations. This stems from the 

A similar plot was presented in iTrac. Bode fc Ostrikei] feoill) , 
except that here the variable is the maximum mass, not the min- 
imum mass. 
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Figure 7. The value of the integral in equation (0 for various 
Zmax (x-axis) and M max (y-axis) relative to the fiducial values 
Zrnax = 2 and Mmax = 10 16 h~ 1 M Q , for I = 4000. Contours 
show constant fractions of the fiducial value; the calculations were 
performed using the fitting formula in equation 114> . 

appreciable impact of baryonic processes (both stellar and 
gaseous) on cluster evolutio n. As an example we briefly men - 
tion here the conclusion of IStanek. Rudd fc Evrardl (|2009| ) 
that number densities of clusters in their hydrodynamical 
simulations deviate by 10% — 60% from those predicted by 
iTinker et all i|2008r i. 



4 DISCUSSION 

We have presented a new analytical method of calculating 
the thermal SZ power spectrum and cluster number counts. 
This merger-tree method takes into account the intrinsic 
scatter in the halo parameters and their dependence on the 
redshifts of formation and observation that results from the 
hierarchical evolution of clusters. In particular, our model 
provides a good statistical description of the high-mass and 
high-redshift halos. 

We have demonstrated that our approach, which does 
not rely on pre-calibrated scaling relations, allows to explore 
different density profiles of DM halos. In addition, since the 
cosmological parameters can be changed in each run of the 
merger-tree code and no scaling relations are inserted by 
hand, our approach is more suitable for cosmological param- 
eter estimation studies than the standard analytic method. 

We stress that the approach presented here differs from 
adding scatter in the various observable parameters to the 
standard analytic calculation, since in this case the scaling 
relations and the scatter itself are inserted by hand, whereas 
here they are intrinsic to the model and its predictions. 

Our calculations predict a higher power spectrum nor- 
malization than the standard ca l culatio n which uses the 
c(M, z) relation from iDuffv et all (|2008h (for the same IC 
gas model), and a shift of the peak to higher multipoles. In 
addition, the merger-tree model predicts higher cluster num- 
ber counts than the standard approach (again, for the same 
IC gas model). These differences can be explained by the 
modified c(M, z) relation that is predicted by the merger- 
tree model. In particular, it has a higher normaliz ation and a 
slower evolution with redshift than the results of lDuffv et alj 



(2008). We argue that the higher normalization of the con- 
centration parameter is more compat ible with ob s erva tions, 
while the faster redshift evolution of lDuffv et al.l l|2008h re- 
sults from the fact that their sample is heavily weighted by 
low-mass low-redshift halos, where the dependence on red- 
shift is expected to be stronger. The excellent agreement 
between ou r results and the adiab atic hydrodynamical sim- 
ulations bv lBattaglia et all (|2010h further support our con- 
clusions. 

As anticipated, we have found the SZ power spectrum 
is quite sensitive to the assumed IC gas model. Among the 
parameters that can strongly influence the SZ power are the 
gas mass fraction and the amount of non-thermal pressure 
and their variation with cluster mass and redshift. Our fidu- 
cial model, which assumes a constant f g = 0.1 produces a 
level of the SZ power spectrum significantly higher than re- 
cent measurements by the SPT. When we adopt a scaling 
of the gas ma s s frac tion with mass f g oc M - 2 , inferred by 
iGiodini et alj (120091 ') . the (current) SPT measurement falls 
within the range of uncertainty in as which extends to lower 
power levels due to the mass-dependent gas fraction. Clearly, 
more detailed modeling of the gas mass fraction is needed 
for meaningful comparison with observations of the thermal 
SZ effect. Note also that we adopted a class of poly tropic 
gas models; different models can be readily explored in the 
context of our merger-tree approach. 

Another basic source of uncertainty is the halo MF. 
We have shown that the high-mass, high-redshift halos, for 
which the MF calibration is least certain, have a relatively 
large contribution to the SZ power spectrum. In this work 
we have used the Sheth-Tormen MF, but our approach can 
be extended to include other mass functions. This can be 
achieved by either calibrating the conditional MF that we 
use to build the merger t rees to an up-to-date numerica l 
simulation, in the spirit of iParkinson. Cole fc Hellvl {2008), 
or constructing a theoretically motivated conditional MF 
that provides a better fit to simulations. 

In the calculations presented here we assumed spherical 
symmetry for all halos. However, it is possible to extend our 
approach to account for halo triaxiality by modeling each 
halo in the tree as an ellipsoid either with constant values 
of the axis ratios, or with values drawn randomly from a 
probability distribution. 

In this work we calculated the power spectrum of the 
SZ thermal component. Using our approach it is possible to 
calculate also the power due to the SZ kinematic component 
due to random cluster motions by assigning a peculiar veloc- 
ity to each halo. This contribution is expected to be a small 
fraction of that from the thermal component. However, we 
cannot use our method to calculate the full kinematic power 
which also has contributions from matter outside of clusters. 
Therefore, our approach is better suited for calculations of 
cluster number counts and cosmological parameter estima- 
tion which relies on cluster counts, since they depend only 
on the signal that originates from matter within clusters. 

The approach presented in this work provides a com- 
putationally efficient way to explore the uncertainties dis- 
cussed above, and to gain an improved understanding of the 
influence of the evolution of galaxy clusters on their various 
observable properties. 
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